Thermo-Fluid Modeling of the Pressurization and 
Draining of a 1000-Gallon Cryogenic Tank with 
GFSSP 
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GFSSP 


GFSSP is the Generalized Fluid System Simulation Program. 


It is a general-purpose computer program to compute 
pressures, temperatures, and flow rates in a flow network. 


It was primarily developed to analyze: 
— Internal Flow Analysis of a Turbopump 
— Transient Flow Analysis of a Propulsion System 


GFSSP has been in continuous development at MSFC since 
1994. 


GFSSP 


¢ GFSSP discretizes a system into a flow network of nodes 
connected by branches. 
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GFSSP 


Conservation of mass and energy is solved in the nodes to get 
pressures and temperatures. 


The momentum equation is solved in the branches to get flow 
rates. The branches represent empirical laws of transport 
processes, such as pressure drop in a pipe. 

Integrated property packages (GASP/WASP and GASPAK) 
provide fluid properties. 

Built-in options include 

— Pressure and flow regulators 
Opening/closing of valves 
Heat exchangers 
— Tank pressurization 


Tank Pressurization Option 


¢ GFSSP has a built-in tank pressurization option. 


¢ Asingle ullage node is separated from a propellant node by a 
pseudo-boundary node. 
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Tank Pressurization Option 


e As the model runs and the tank drains, the 





code automatically tracks: ae 
v 
— Volume of ullage and propellant nodes || 
— Propellant depth and head (pgh) as cas 
— Ullage-to-wall heat transfer 5. af 
— Ullage-to-propellant heat transfer \ me 
— Average temperature of uncovered tank wall a ne 


— Pressure and temperature of ullage node 


¢ Heat transfer calculations are based on 
natural convection correlations 


recommended by Elliot Ring (1964). Ke 
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¢ Correction factors may be applied to the : 
heat transfer rate. Ve 


Propellant to Engine 


Tank Pressurization Option 


An early application was modeling of LOx tank pressurization 
during FASTRAC engine testing at Stennis Space Center. 
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IVF 1000 Tank 


e Integrated Vehicle Fluids (IVF) is a test 
program under NASA’s eCryo project, which 
evolved from the Cryo Propellant Storage 
and Transfer (CPST) project. 


¢ The IVF-1000 test tank was formerly known 
as the CPST Engineering Design Unit (EDU). 
— Dimensions: 87” x 67” 
— Volume: 144.5 ft? (1080 gal) 
— Mass: 400 |b of aluminum alloy 
— Insulation: 1.25” of SOFI + 60 layers MLI 
¢ In September 2015, a series of pressurization 
tests were performed on the IVF-1000 tank. 
The tank was filled with LH, and LN,. 





Model Details 


Helium Pressurization System A Bang-Bang Valve 
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Ss Drain Valve 


Helium press system includes bang-bang valve and control orifice. 


Propellant drain valve is closed during pre-press. During drain, it opens to an area chosen to 
match the observed drain rate. 


LH, Test Description 


LH2 propellant 

He pressurant 

Bang-bang valve set to maintain 31-32 psia range. 
Tank drains from 89.2% to 69.8%. 

Pre-press starts at ~T-392 sec. 

Drain starts at ~T+7 sec. 

Drain stops at ~T+295 sec. 
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Tank Fill Level by Capacitance Probe 


—< CP2440 Percent Capacitance Probe 


92 WinPlot v4.55.1 


88 


84 





80 
Initial: 89.2% @ 7 sec 
Final: 69.8% @ 295 sec 


76 
Estimated drain flow rate: =19.4% * (144.5 ft3) / 288 s 
= 0.0973 ft3/s = 43.7 gpm 
72 
= 0.434 Ib/s 
68 
0 50 100 150 200 250 300 


Time Seconds Relative Time in Seconds 


12:14:37PM 12/14/2015 


Results of 15° Run 


Ullage pressure rises more quickly in model. 
Ullage pressure overshoot is much less in model. 
Decay of ullage pressure is much faster in model. 





Therefore, many more pressure cycles in model. 


Difficult to compare GFSSP’s temperature in the single ullage node to the multiple 
ullage temperature measurements, but model ullage temperature appears to be 
“in range”. 

Model pressure upstream of orifice is “10% higher than test measurement. 
Suggests model flow rate is ~10% higher. (Some uncertainty because orifice C, is 
unknown. Assumed 0.84.) 

Model pressure downstream of orifice decays much more quickly than measured 
pressure. Suggests gas in downstream press lines takes its time reaching the tank. 
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Pressure at Orifice 
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Note: Pressure decay is much slower in test 


Run #2 


¢ Aseries of parametric runs were carried out to investigate the effect of various 
model inputs on ullage pressure rise and decay. 


¢ Changes: 
— Control valve C, reduced from 0.6 to 0.5. 
— Set heat transfer correction factor to 0.1. 
— Add 200 in? of extra volume to represent press line to aft diffuser. 
— Set temperature of pressurant to match measured diffuser temperature. 


Init. P,,, (psia) | Init. Press Time (s) | Init. a (CELA) a 


Run #1 33.9 
Run #2 32.8 13.5 75 32 
Test Data 38.5 5.4 140 24 
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Pressure at Orifice 
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LN, Test Description 


LN2 propellant 


He pressurant 


Bang-bang valve set to maintain 37.5 — 38.5 psia range during pre-press, 
and 31 — 32 psia range during drain. 


Tank drains from 90.3% to 69.9%. 
Pre-press starts at ~T-1032 sec. 


Drain starts at ~T-O sec. 
Drain stops at ~T+864 sec. 


Adjustments developed for LH2 model are retained in LN2 model: 


Control valve C, reduced from 0.6 to 0.5. 
Set heat transfer correction factor to 0.1. 
Add 200 in? of extra volume to represent press line to aft diffuser. 


Set temperature of pressurant to match measured diffuser 
temperature. 
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Tank Fill Level by Capacitance Probe 
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Initial: 90.3% @ 0 sec 
Final: 69.9% @ 864 sec 
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Results of 1°* Run 


Rate of ullage pressure rise in model is similar to test. 

Ullage pressure overshoot is much less in model. 

Decay of initial ullage pressure (when tank isn’t draining) is much slower in model. 
But after a few cycles, decay of ullage pressure is faster in model than in test. 


Difficult to compare GFSSP’s temperature in the single ullage node to the multiple 
ullage temperature measurements, but model ullage temperature appears to be 
“in range”. 


Model pressure downstream of orifice decays much more quickly than measured 
pressure. Suggests gas in downstream press lines takes its time reaching the tank. 
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Note: Pressure decay is much slower in test 


Discussion 


¢ Heat transfer is the greatest uncertainty in pressurization modeling. 


— The natural convection correlations recommended by Ring (1964) significantly 
over-predicted the heat transfer for this tank and diffuser. 


— Past modeling of FASTRAC tests did not require this correction factor. Perhaps 
the large Stennis facility tank better approximates a flat plate. 

— The single-node ullage approach is an approximation which doesn’t capture 
temperature stratification in the ullage. Results may be better with a well- 
mixed ullage. 

¢ Uncertainty in the total volume of the press line also appears to be significant. 
— Test data always shows more overshoot than predicted by model. 


— Test data shows facility line pressure decays much more gradually than 
predicted by model. 


e There were no measurements of pressurant temperature just as it enters the tank. 
This would be desirable in future testing. 


¢ Almost all of the pressure drop in the press system occurred in the valves and 
control orifice. Pipe roughness and minor losses from bends had little effect. The 
press system portion of the model could probably be greatly simplified. 


Forward Work 


The project has continued with pressurization testing of a larger (4500 gal) tank, 
with a flight-like diffuser (spring 2016) and an injector (Summer 2016). 


We are exploring multi-node ullage set-ups to approximate stratification and 
mixing: 
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We will explore whether the single-node ullage heat transfer is better modeled 
when there is strong mixing induced by an injector. 


